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Resumo 

Nesse pequeno trabalho discutimos estados de equilibrio de estrealas, consideradas como um corpo de 
gas sujeito a agao da forga gravitacional propria, utilizando um modelo analitico simplificado. Usamos uma 
condigao para o estado de equilibrio de um corpo de gas na forma de uma equagao diferencial que relaciona 
distribuigao de pressao e densidade de massa no corpo. Foram discutidas as distribuigoes de densidade de 
massa da forma constante, potencial, exponencial e gaussiana. Foram obtidas as expressoes exatas para 
distribuigao de massa e pressao na diregao radial e a pressao central. 

Palavras-chave: modelo analitico estrelar, estado de equilibrio, solugoes exatas. 

Self-graviting Gas Spheres in Equilibrium State 

Abstract 

In the paper we discuss equilibrium states of stars, using a simplified analytic model. A star is considered 
as self-graviting body of gas. We use a condition for the equilibrium state of the body in the form of a 
differential equation, which relates the pressure distribution and mass density in the body. The density 
distributions of constant, potential, gaussian, and exponential forms are discussed. Exact expressions for 
the distribution of mass and pressure along the radial direction, and the central pressure were obtained. 

Key-words: stellar analytic model, equilibrium state, exact solutions. 



1 Introdugao 

De acordo com uma visao moderna uma estrela e considerada como uma esfera massiva de gas (ou 
mais exato de plasma), mantida Integra pela forga de gravitagao propria. A estrutura interna de 
estrelas e bastante complicada e depende de massa estrelar. Para estrelas de massa solar (0,5-1,5 de 
massa solar) da sequencia principal da diagrama de Hertzsprung-Russell a estrutura interna inclui; 
um nucleo, uma zona de radiagao e uma zona de covecgao. As estrelas massivas (com massa maior 
de 1,5 de massa de Sol) possuem nucleo convectivo acima de qual e localizada a zona de radiagao (ou 
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um envelope de radiagao). Nas estrelas de massa menor de 0,5 de massa de Sol a zona de radiagao e 
ausente (diz-se que estrelas possuem nm envelope de convecgao). Nas estrelas da seqnencia principal 
a energia esta sendo gerada pela queima de hidrogenio em helio atraves de fusao nuclear em sen 
niicleo. 

Na classificagao de estrelas fora da seqnencia principal sao distinguidas estrelas de tipos: sub- 
gigantes, gigantes, gigantes luminosas, supergigantes, liipergigantes, sub-anas, estrelas retardatarias 
azuis. Nas estrelas destes tipos a energia e gerada principalmente pelo mecanismo de fusao nuclear. 
Alem disso sao consideradas as estrelas prAsequencia principal. A fonte de energia desses objetos e 
causada somente pela contragao gravitacional em oposigao a fusao nuclear em estrelas de outros tipos. 
Estrelas de varios tipos possuem caracteristicas especfficas de estrutura, mecanismos da produgao e 
transporte de energia. Para alguns tipos de estrelas sao construidos os bons modelos que descrevem 
adequadamente as caracteristicas observadas das estrelas, para outros tipos os modelos estao em 
processo de construgao e discussao. 

O modelo mais simples de estrutura estelar e a aproximagao quase-estatica de simetria esferica. 
Uma discussao desse modelo e apresentada em Refs. n. 0. O modelo assume que a estrela possui 
simetria esferica e esta em estado de equilibrio hidrostatico. O modelo e baseado em um conjunto de 
equagoes diferenciais ordinarias. Duas equagoes descrevem variagao de materia e pressao na diregao 
radial. Outras equagoes descrevem variagao da luminosidade e transporte de energia. Para deter- 
minagao da luminosidade, precisa da caracteristica taxa de produgao de energia. Esta caracterfstica 
e determinada atraves de um mecanismo principal de transporte de energia: radiagao on convecgao. 
O conjunto das equagoes do modelo compoe um sistema das equagoes nao lineares, que pode ser 
resolvido em geral por metodos numericos. 

De fato nesse modelo as duas equagoes para materia e pressao sao desacopladas das outras 
equagoes, pois elas nao contem parametros como temperatura, taxa de produgao de energia, lu¬ 
minosidade, portanto podem ser consideradas separadamente. Impondo uma fungao da distrubuigao 
de densidade da materia na estrela, e possivel obter a distribuigao de pressao. Adicionando uma 
equagao de estado para parametros termodinamicos, como por exemplo a equagao de estado de um 
gas perfeito, e possivel calcular a temperatura na estrela. Adicionando ainda uma relagao apropriada 
para taxa de produgao de energia, e possivel estimar a luminosidade da estrela. De tal maneira pode 
ser construido um modelo analitico. Claro que esse modelo e simplificado, mas permite fazer es- 
timagoes de algumas caracteristicas principais de estrela, tais como: densidade central, temperatura 
central, luminosidade em termos de raio e massa de estrela. Tal forma de construgao de modelo 


2 


analitico foi discutida em Ref. [3], onde foi usada a distribuigao de densidade na forma linear, cha- 
mado modelo linear de estrela. Entretanto, o nso da densidade na forma linear nao parece mnito 
realistico. 

Nesse peqneno trabalho nsamos os modelos analiticos com varias distribnigoes de densidade com 
aspectos mais realisticos. Uma pecnliaridade dos modelos propostos no trabalho e a possibilidade de 
resolver as eqnagoes diferenciais analiticamente. No trabalho focalizamos na obtengao da distribnigao 
de pressao na diregao radial e pressao central, ignorando a determinagao de ontros parametros. Dis- 
cntimos os modelos com distribnigoes de densidade das formas potencial, exponencial e ganssiana. No 
inicio discntimos tambem nm modelo com densidade constante para demonstrar transparentemente 
os passos realizados nos calculos, e para introduzir notagoes nsadas posteriormente. As distribnigoes 
de massa e de pressao obtemos na forma anahtica. Para realizar os calcnlos foram usados varios 
metodos de fisica matematica disponiveis na literatnra, por exemplo 0, 0. Para visnalizagao e 
comparagao dos resnltados apresentamos tambem os resnltados grahcamente. Na conclusao discnti¬ 
mos os resnltados obtidos e snas aplicagoes possiveis. 

Um dos aspectos motivacionais deste trabalho e a sna aplicagao pedagogica. No trabalho demos- 
tramos de forma concisa a constrngao das eqnagoes e da condigao de contorno, qne descrevem nm 
modelo analitico e simplihcado das estrelas. Isto e, a eqnagao de eqnilibrio hidrostatico e a eqnagao 
de conservagao de massa. Snpomos qne os modelos discntidos neste trabalho podem ser ntilizados 
em salas de anla on como exercicios em disciplinas introdntorias de astrofisica, disciplinas dedicadas 
a modelagem de fenomenos fisicos, on ate mesmo em eqnagoes diferenciais ordinarias. 

2 Esferas de Equilibrio 

Consideraremos nm volnme de gas esfericamente simetrico de raio R em estado de eqnlibrio. A 
coordenada radial r dentro da esfera varia como 0 < r < R. Escolhemos nm fragmento de volnme 
dentro da esfera da forma de nm paralelepipedo retangulo inhnitesimal dV = dr As com base As e de 
altura dr (ver Fig. [T]). No fragmento de volnme atnam as segnintes forgas externas: a forga snperhcial 
de pressao das camadas adjacentes de gas e a forga gravitacional volnmetrica. Para garantir eqnilibrio 
do fragmento dentro da esfera de gas a forga gravitacional Fg deve ser eqnilibrada pelas forgas de 
pressao qne atnam na face inferior do fragmento Fpi e na face snperior do fragmento Fp 2 (forgas 
de pressao perperndicnlar a diregao radial sao eqnilibradas por cansa de simetria esferica e nao sao 
indicadas na hgura), portanto: 

Fg + Fpi -|- Fp2 = 0 , ( 1 ) 
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onde: 


^ ^Mrdm „ , „ 

Fg = G - Fpi = Pi As, Fp 2 = P 2 AS , 


( 2 ) 


dm = pdV = pdrAs e a massa do fragmento infinitesimal considerado, pea densidade da substancia 
de gas, Pi = p (r), p 2 = p (r + dr) sao pressoes nos pontos indicados, Mr e a massa dentro da esfera 
de raio r 


Mr = / 47rp (r) r'^dr 

Jo 


Levando em conta os sentidos das forgas, temos: 

^MrpdrAs 


G- 


+ p (r) As — p (r + dr) As = 0 . 


(3) 


(4) 


Apresentando a diferenga das pressoes nos pontos proximos: 

p (r) — p (r + dr) = — ^dr, 

dr 


(5) 


chegamos a equagao diferencial ordinaria 


^ 

dr 


( 6 ) 


que relaciona a fungao de distribuigao de densidade da substancia de gas p (r) e pressao p (r). A Eq. 
(E]) apresenta a equagao de Euler para um fiuido em estado de equilibrio para o modelo considerado. 
Na superficie da esfera a pressao e nula. Obtendo assim a condigao de contorno no ponto r = R: 


p (R) = 0 . 


(7) 


Mostraremos varias formas de distrubuigao de densidade da substancia na esfera que faz com que 
obtenhamos solugoes exatas para distribuigao da pressao. 


2.1 Desidade Constante 

Escolhemos que a densidade da substancia dentro da esfera seja constante 

p (r) = po = const . 

Pela Eq. ([3]) temos para Mr'. 

Mr = I dvrpor^dr = ■^por'^ = 


" 2^. _47r_3 AnpoR^ fr\5 


3 \R 


sendo que a massa total da esfera e: 


Mr = Mr {R) = ^R^ . 


( 8 ) 


(9) 


( 10 ) 
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Figura 1: Esquema das forgas externas atuando no fragmento de volume infinitesimal. 



Para determinagao de p (r) da Eq. ([6]) temos a equagao diferencial: 


dp 

dr 






Po 


on 

^ ^ inGpl r 
dr 3 i? ’ 

com a condigao de contorno ([7]). Resolvendo a Eq. fITT]) . temos: 


P (r) -po = 


3 Jo R 



2I 

3 2 \rJ 


( 11 ) 
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P/Po> Mj/Mr 



P/Po 



Figura 2: Distribuigoes de p, Mr, p para p = const. 


onde po = p (0) e pressao no ponto r = 0, que apresenta a pressao cetral na esfera de gas. Aplicando 
a condigao de contorno ([7]), obtemos: 


Po=2 


1 dvrGpo 




( 12 ) 


portanto: 


p (r) = 


AirGpl 1 


^2 \ 


As distribuigoes de p(r), Mr (r), p (r) sao mostrados na Fig. [2j 
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2.2 Fungao de Distribuigao da Densidade na Forma Potencial 

Escolhemos a distribuigao de densidade da substancia na forma potencial 

• ( 14 ) 
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Figura 3: Distribuigoes de p, Mr, p para p na forma potencial com n = 2,5, 8. Nos graficos a curva 
que corresponde a n = 5 e intermeiaria entre as curvas den = 2e?7, = 8. 


Pela Eq. ([3]) temos para Mp. 


Mr= dvrpo 1 - ^ r'^dr = 




dvrpo-R^ 

3 R 


3 M 
n + 3 i?" 


sendo que a massa total da esfera e: 


Mr = Mr {R) = 


dvrpo-R^ n 


3 n + 3 

Para determinagao de p (r) da Eq. ([6]) temos a equagao diferencial: 

dp AttGpqR ^ a n + 6 r"" 3 

dr ~ 3 V n+ 3 R^ ^ n + 3R‘^^ 


Resolvendo a Eq. on, temos: 
p(r) -po = - 


2R2 


2 (n + 6) r" 

1 - ^ ^ — + 


^2n 


{n + 3) (n + 2) R" {n + 3) (n + 1) R^^ 


Aplicando a condigao de contorno ([7]), obtemos: 


”32 


^ 2(n + 6) 


(n + 3) (n + 2) (n + 3)(n + l)J ’ 


(15) 


(16) 


(17) 


(18) 


(19) 
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portanto 


p (r) = 




'2 1 
3 2 


1 - 


2 (n + 6) 


+ 


(n + 3) (n + 2) (n + 3) (n + 1) 




2 (n + 6) R 


+ 


3 


( 20 ) 


{n + 3) (n + 2) R"^ {n + 3) (n + 1) R^'^ 

As distribuigoes de p(r), Mr{r), p{r) sao mostradas na Fig. |3] para n = 2, 5, 8. No grafico as 
distribuigoes de (r) sao fungoes crescentes. Observamos tambem que na fungao de distribuigao 
de densidade d a potencia n nao e necessariamente um niimero inteiro e pode ser qualquer real 
n > 1. 


2.3 Fungao de Distribuigao da Densidade na Forma Exponencial 


P/Po> Mj/Mr 



P/Po 



Figura 4: Distribuigoes de p, ilA-, p para p na forma exponencial. 











Escolhemos a distrubuigao de densidade da substancia na forma: 


P [r) = po 


1 - 


3R 


e « , 


sendo que 


Pela Eq. ([3]) temos para 


P (R) = Po 


1 

1 - - 
3 


Mr = dvrpo 

= I 


r 

1 - — 

Jo 

L 3R\ 


'0 


1 - — 
3R. 


1 2 

- — irPo 
e 3e 


e Rr^dr 




uo 




= AtipqR^ 

Na Eq. fl25D . usando a tecnica de integragao por partes: 

^ \ 3 


'0 




R 


e-,+3l . 


ry \ Z / T’ 


portanto: 




sendo que a massa total da esfera e: 


Mr = Mr (R) = ^R^ , 
3e 


que permite tambem escrever Mr como: 


^ \ 3 


Mr = MR{-] e^-R . 




RJ 


Para determinagao de p (r) da Eq. ([6]) temos a equagao diferencial: 

dp AirGplR ^ A 1 r \ _2j_ 

dr 3 R\ 3 r) 


Resolvendo a Eq. fl28|) . temos: 

P (r) -Po = 
AttGpIR^ 


^nGplR^ 

3 .L R\ 3 R 


r \ _r,r_ ,/r\ 1 r / r \ ‘^ 


L-'o 


-]e-‘^^d(^ 


RJ 3 ./o VR/ 


Na Eq. (j29]): 


portanto: 


r-r , „ . 2 _ / ^ 


.-2^. 


R 


1 /r \2 


2 \R 


AttGpIR^ 


P {r) -po = 


2 r 
sio 


— 2^ 
o ^ D 




( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


\R 




^ r \' 

VR/ 

(29) 


(^) 

. (30) 

r \ 2 

—2^ 
e ^ 

(31) 
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A integral 




r 1 

^^2/^4’ 


portanto 


P (r) -Po = 


AttGpIR^ 


6 \ 6 3R 6\rJ 


Aplicando a condigao de contorno ([7]), obtemos: 

,47rpo-R^ /I 1 


Po = G- 


6 3e2 


portanto 


AnGplR^ 


p [r) = 


1 1\ /I / I Ir l/r\ 2,2 

6“^j “ V6^ V~6 ^ 


47rG'p^i?^ 1 


1 + ^-1 
2 R 2 \rJ j 


3 3e2 

As distribuigoes de p(r), (r), p{r) sao mostrados na Fig. |H 


(32) 


(33) 


(34) 


(35) 


2.4 Fungao de Distribuigao da Densidade na Forma Gaussiana 

Escolhemos a distrubuigao de densidade da substancia na forma gaussiana: 

2 


sendo que: 


Pela Eq. ([3]) temos para Mp. 


P [r) = Po 


P i.R) = Po 


1 - 


3i?2 


e , 


2 

1 - - 
3 


Mr = dvrpo 
= dvrpo^^ [ 


pr 

r 2 r 2 1 
1 

L 

3i?2 


'0 L 


2 


1 1 

“ — irPo ■ 
e de 


e T^r^dr 


^ Rii 


= dvrpo-R^ 


uo 




Observando que: 


de ^ = —2xe ^ dx , 


temos para a integral: 


r \ 4 / r 

' " \r) [r 


1 

2 


7 > N 3 _ ^2 

-) * 7!. = 


3 _ ^2 

e 




0 -/o 




r\3 r _.2 

— )e«2^—3 / 

R' Jr\ 


r \ 2 / r 


(36) 


(37) 


(38) 


(39) 


( 40 ) 
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P/Po. Mj/Mr 



P/Po 



Figura 5: Distribuigoes de p, Mr, p para p na forma gaussiana. 


que da: 


Mr 


A^pqR^ 







(41) 


Portanto: 


Mr 


47rpo-R^ 

3 



sendo que a massa total da esfera e: 


Mr = Mr {R) = ^R^ , 

6 e 

que permite tambem escrever Mr como: 

/ r \ 3 , r2 

Mr = Mr e^-1^ . 


(42) 


(43) 


(44) 
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Para determinagao de p (r) da Eq. ([6]) temos a equagao diferencial: 


Resolvendo a Eq. temos: 


dp 

AttGpI r 

1 - 

2 r 2 ' 

e- 2 - 

dr 

3 R 


3 R 2 J 



AnGplR^ 1 



2 r 2 

Po - 

- 3 i 

R 

3R\ 


e ^ d— 

-/l 


( 45 ) 


r r 


2 rr ^3 


—e d— — - I -^e "^^d— 


uo 


R 


R 3 /n R3 


R 


(46) 


Na Eq. fl46l) a integral: 


f'’ .^3 
' -2 


R3 


r2 r 1 


r , 

= — 



R A 

.J 



1 


—-e + 



R2 

2 

\ 


- 2 ^ 


1 

'4 

1 

“4 


7^2 2 

^ - 22 :^ 

e ^ 


R 2 


2 7^2 

- 22 ^ 7 

e ^ d—— 


R 2 


/ 1 \ 2 ''^ 1 
-h - ^- 

\R^ 2) 2 


e a integral 


(47) 


Jo di 
Portanto: 


r .24, r 1 r -24jr^\ 1 f -24,/-2r^ 

—e ^d — = - e ^d —^ = — / e ^d 


R 
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\R^ J 


'0 


R 2 


= -- ( - 1 


(48) 


P (r) -Po = 


AtiGpIR^ 


\ -p 4 


-h - I e“^:^ 

R2 ^ 2 ' 


= -G 


47rp2i?2 1 

3 6 


/ ^2 




Aplicando a condigao de contorno ([7j), obtemos: 


(49) 


Po = G 


Anp^^R^ 1 
3 6 


(50) 


Portanto 


p (r) = 


f 1 _ AP-# 
3 el IP 


As distribnigoes de p (r), (r), p (r) sao mostrados na Fig. [5l 


(51) 


3 Conclusao 

Neste trabalho foram discntidos modelos estrelares analiticos simplificados que inclnem dnas eqnagoes 
basicas: a eqnagao de eqnilibrio liidrostatico e a eqnagao de conservagao de massa. Usando distri- 
bnigao de densidade da massa na estrela como nma determinada fnngao obtemos analiticamente 
a correspondente distribnigao de pressao dentro da esfera e determinamos a pressao central. Dis- 
cntimos as distribnigoes de densidade da forma constante, potencial, exponencial e ganssiana. As 
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distribuigoes propostas no trabalho sao ate um determinado gran convencionais. Entretanto elas sao 
mais realisticas (exceto o exemplo de densidade constante) comparando, por exemplo, com o modelo 
estrelar linear. 

Para modelagem de estrelas com nm micleo extenso e mais apropriada a distribnigao da forma 
potencial; Variagao do gran da fnngao potencial possibilita modelar a extensao do nncleo estrelar. 
Para modelagem de alguns tipos de estrelas e mais apropriado ntilizar nm modelo com densidade 
na snperficie nao nnla. As distribnigoes exponencial e ganssiana realizam esse mesmo cenario. Em 
algnns casos pode ser conveniente fazer estimagoes, utilizando a condigao de contorno no infinito. 
Os modelos com as distribnigoes exponencial e ganssiana podem ser facilmente reformnlados para a 
condigao de contorno no infinito: p (r) = 0 qnando r —> oo. No caso da distribnigao exponencial a 
correspondente solngao toma a forma: 


p (r) 


AttGpIR^ 1 
3 6 




(52) 


com a pressao central 


Po = 


At^Gp^R^ 1 


3 6 

No caso da distribnigao ganssiana a solngao tern a mesma forma das Eqs. (Ea, dsiii 


(63) 


p (r ) = 


AtiGpIR^ 1 
3 6 


\ 


-2 




1 Po — G 


3 6 


(54) 


Para dar nm exemplo de nma distribnigao da densidade dentro de estrelas reals apresentamos 
na figura ([S]) a distribnigao de densidade dentro do Sol obtida pelo Birmingham Solar Oscillations 
Network (BiSON-13) da Ref. [7]. Prestamos atengao qne os dados na Ref. iniciam do valor 
de t/Rq = 0,058884 {Rq e o raio do Sol). Portanto para constrnir o modelo de qnalqner forma 
e necessario aplicar nma hipotese sobre o intervalo central 0 < r < Para estrelas distantes e 
impossiVel obter a informagao tanta detalhada. Usando os dados observacionais diretos e indiretos, 
fazem conclusoes sobre as caracteristicas mais gerais: Inminosidade, massa, raio, caracterfsticas 
espectrais, etc.. Entao, modelagem fisica matematica torna-se o linico metodo de investigagao da 
estrntnra interna das estrelas. 

As fnngoes propostas neste trabalho a titnlo de distribnigoes da densidade podem ser consideradas 
como fnngoes de teste. Introdngao de um parametro (on varios parametros) na fnngao possibilita 
ajustar a fnngao a um aspecto mais apropriado pela variagao do parametro. Se ainda satisfazer a 
propriedade de obtengao de solugoes das equagoes do modelo na forma anahtica podemos chegar a um 
modelo (de fato a um classe de modelos) analitico exatamante soluvel com possibilidade de escolha 
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Figura 6: Distribugao da densidade de massa no Sol segundo os dados de BiSON-13 [7]; a densidade 
p e apresentada em g/cm^. 

e ajuste da fungao de distribuigao da densidade. Tanta discussao de modelos analiticos exatamante 
soliiveis vemos como uma das possiVeis continuagoes deste trabalho. 
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